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We consider pion electroproduction on a proton target close to threshold for 
Q 2 in the region 1 — 10 GeV 2 . The momentum transfer dependence of the 
S-wave multipoles at threshold, -Bo+ and Lq+i is calculated using light-cone 
sum rules. 
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1. Threshold Pion Production 

Pion electroproduction at threshold from a proton target 

e(i) +p(P) -» e(l') + 7T+(k) + n(P') , 
e(l)+p(P)^e(l')+n (k)+p(P') (1) 

can be described in terms of two generalised form factors defined as 1 

(N(P')n(k)\j; m (0)\p(P)) = (2) 

= -f iV(P')75 {M q, 4) — 2 Gl N m - i -^f G r(Q 2 )) N(P) , 

Jit ^ ZTYl J 

which can be related to the S-wave transverse Eq + and longitudinal Lq+ 
multipoles: 



N _VWa^ I (2m + m^) 2 + Q 2 f 2 N 1 N 



N _ V47ra cm mm h | / (2m + m^) 2 + Q 2 ( N 2m n N \ 
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The differential cross section at threshold is given by 



da. 



7* 



2\k f \W 



th W 2 



(^) 2 + ^«) 2 • (4) 



Here and below m — 939 MeV is the nucleon mass, W 2 = (k + P 1 ) 2 is 
the invariant energy, kf and are the pion three-momentum and the 
photon energy in the cm. frame. The generalised form factors in (2) are 
real functions of the momentum transfer Q 2 at the threshold W = m + . 
For generic W the definition in (2) can be extended to specify two of the 
existing six invariant amplitudes, Gi : 2(Q 2 ) — > Gi.2{Q 2 , W), which become 
complex functions. 

The celebrated low-energy theorem (LET) 2-4 relates the S-wave multi- 
poles or, equivalently, the form factors G\, Gi at threshold, to the nucleon 
electromagnetic and axial form factors for vanishing pion mass = 

Q 2 ^7r°p _ 9 A Q 2 nP nn ° p _ 2g A m 2 np 

m 2 ^ 1 ~2(Q 2 + 2m 2 r M ' ^ 2 ~ (Q 2 + 2m 2 ) ' [ ) 

_ \/5 (Q 2 + 2m 2 ) \71 ' 2 " (Q 2 + 2m 2 ) GE - 

Here the terms in Gm,e are due to pion emission off the initial proton state, 
whereas for charged pion in addition there is a contribution corresponding 
to the chiral rotation of the electromagnetic current. 

The subsequent discussion concentrated mainly on the corrections to (5) 
due to finite pion mass 5,6 . More recently, the threshold pion production for 
small Q 2 was reconsidered and the low-energy theorems re-derived in the 
framework of the chiral perturbation theory (CHPT), see 7 for a review. 
The new insight gained from CHPT calculations 8 is that the expansion at 
small Q 2 has to be done with care as the limits — > and Q 2 — > do not 
commute, in general. The LET predictions seem to be in good agreement 
with experimental data on pion photoproduction 9 , However, it appears 
10,7 that the S-wave electroproduction cross section (4) for already Q 2 ~ 
0.1 GeV 2 cannot be explained without taking into account chiral loops. 

For larger momentum transfers the situation is much less studied as 
the power counting of CHPT cannot be applied. The traditional derivation 
of LET using PCAC and current algebra does not seem to be affected as 
long as the emitted pion is 'soft' with respect to the initial and final state 
nucleons simultaneously. The corresponding condition is, parametrically, 
Q 2 <C A 3 /m,r (see, e.g. 5 ) where A is some hadronic scale, and might be 
satisfied for Q 2 <~ 1 GeV 2 or even higher. We are not aware of any dedicated 
analysis of the threshold production in the Q 2 ~ 1 GeV 2 region, however. 
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It was suggested in Ref. 11 that in the opposite limit of very large mo- 
mentum transfers the standard pQCD collinear factorisation approach 12,13 
becomes applicable and the helicity-conserving G1 N form factor can be 
calculated for = in terms of chirally rotated nucleon distribution 
amplitudes. In practice one expects that the onset of the pQCD regime 
is postponed to very large momentum transfers because the factorisablc 
contribution involves a small factor al(Q)/n 2 and has to win over nonper- 
turbative "soft" contributions that are suppressed by an extra power of Q 2 
but do not involve small coefficients. 

The purpose of this study is to suggest a realistic QCD-motivated model 
for the Q 2 dependence of the form factors alias S-wave multipoles at 
threshold in the region Q 2 <~ 1 — 10 GeV 2 that can be accessible in current 
and future experiments in Jefferson Laboratory and elsewhere (HERMES, 
MAMI). 

2. Light-Cone-Sum Rules 

In Rcf. 14 we have developed a technique to calculate baryon form factors for 
moderately large Q 2 using light-cone sum rules (LCSR) 15 > 16 . This approach 
is attractive because in LCSRs "soft" contributions to the form factors 
are calculated in terms of the same nucleon distribution amplitudes (DAs) 
that enter the pQCD calculation and there is no double counting. Thus, 
the LCSRs provide one with the most direct relation of the hadron form 
factors and distribution amplitudes that is available at present, with no 
other nonperturbative parameters. 

The same technique can be applied to pion electroproduction. In Ref. 1 
the Gi and G2 form factors were estimated in the LCSR approach for 
the range of momentum transfers Q 2 <~ 5 — 10 GcV 2 . For this work, 
we have reanalysed the sum rules derived in 1 taking into account the 
semi-disconnected pion-nucleon contributions in the intermediate state. We 
demonstrate that, with this addition, the applicability of the sum rules can 
be extended to the lower Q 2 region and the LET results in (5) are indeed 
reproduced at Q 2 <~ 1 GcV 2 to the required accuracy 0{m v ). The results 
presented below essentially interpolate between the large-Q 2 limit consid- 
ered in 1 and the standard LET predictions at low momentum transfers. 

Accurate quantitative predictions are difficult for several reasons, e,g, 
because the nucleon distribution amplitudes are poorly known. In order to 
minimise the dependence of various parameters in this work we only use 
the LCSRs to predict certain form factor ratios and then normalise to the 
electromagnetic nucleon form factors as measured in the experiment, see 
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for the details. In particular we use the parametrisation of the proton 
magnetic form factor from 18 and for the neutron magnetic form factor from 
19 . For the proton electric form factor we use the fit 20,18 to the combined 
JLab data in the 0.5 < Q 2 < 5.6 GeV 2 range 

qp 

^7^ = l-0.13(Q 2 -0.04) (6) 

and put the neutron electric form factor to zero, which should be good to 
our accuracy. Note that using (6) for larger values of Q 2 up to 10 GeV 2 is 
an extrapolation which may be not justified. 

The resulting LCSR-based model is shown by the solid curves in Fig. 1, 
where the four partial waves at threshold that are related to the generalised 
form factors through the Eq. (3) are plotted as a function of Q 2 , normalised 
to the dipole formula 

G D (Q 2 ) = 1/(1 + Q7moT (7) 

where /Zq = 0.71 GeV 2 . This model is used in the numerical analysis pre- 
sented below. We expect that its accuracy is about 50%. It can be improved 
in future by the calculation of radiative corrections to the LCSRs, especially 
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if sufficiently accurate lattice calculations of the moments of nucleon dis- 
tribution amplitudes become available. To give a rough idea about possible 
uncertainties, the "pure" LCSR predictions (all form factors and other in- 
put taken from the sum rules) are shown by dashed curves for comparison. 

3. Moving Away From Threshold 

As a simple approximation, we suggest to calculate pion production near 
threshold in terms of the generalised form factors (2) and taking into ac- 
count pion emission from the final state which dominates the P-wave con- 
tribution in the chiral limit (cf. 11 ). In particular, we use the following ex- 
pression: 

(N(P')n(k)\j^(0)\p(P)) = 

= -j/^ {(W i) ^Gl N (Q 2 ) '^f Gl N (Q 2 )} N(P) 

+ 2 ;V+1)V^ (/y) w ' + m) {*? m - ¥) 

+ ^fFi(Q^N(P). (8) 

Here Ff(Q 2 ) and F^iQ 2 ) are the Dirac and Pauli electromagnetic form 
factors of the proton, c^o = 1 and — \[2 is the isospin factor, qa = 1.267 
and ,U = 93 MeV. 

The separation of the generalised form factor contribution and the final 
state emission in (8) can be justified in the chiral limit — > but involves 
ambiguities in contributions <~ 0{m^). We have chosen not to include the 
term ~ }t in the numerator of the proton propagator in the third line in (8) 
so that this contribution strictly vanishes at the threshold. In addition, we 
found it convenient to include the term ~ qp.il q 2 in the Lorentz structure 
that accompanies the F\ form factor in order to make the amplitude for- 
mally gauge invariant. To avoid misunderstanding, note that our expression 
is not suitable for making a transition to the photoproduction limit Q 2 = 
in which case, e.g. pion radiation from the initial state has to be taken in 
the same approximation to maintain gauge invariance. 

The virtual photon cross section can be written as a sum of terms 



(9) 
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with 



\M- r \ 2 =M T + eM L + V2e(l + e) M LT cos(^) 

+ eM TT cos(2^) + \y/2e(l - e) M' LT sin(0 7r ) ; (10) 



in the last term A is the beam helicity. 

The complete expressions for the invariant functions are rather cum- 
bersome but are simplified significantly in the chiral limit — ► and 
assuming kf = ©(m^). We obtain 



spin asymmetry contribution ~ M' LT involves imaginary parts of the gen- 
eralised form factors that arise because of the final state interaction. In our 
approximation Mtt — which is because we do not take into account the 
D- and higher partial waves. Consequently, the ~ cos(2^) contribution to 
the cross section is absent. 

We find that the integrated cross sections scale like t7 , 7 » p _ > ^jv ~ 1 / Q 6 , 
which is in agreement with the structure function measurements in the 
threshold region by E136 21 . The S-wave contribution appears to be larger 
than P-wave up to W ~ 1.16 GcV. The ratio of 7r°p and n + n final states 
is approximately 1 : 2 and almost Q 2 -independent. 

The comparison of our calculation for the structure function F^iW, Q 2 ) 
in the threshold region W 2 < 1.4 GcV 2 to the SLAC E136 data 21 at the 
average value Q 2 = 7.14 GcV 2 and Q 2 = 9.43 GcV 2 is shown in Fig. 3. The 
predictions are generally somewhat below these data (~ 50%), apart from 
the last data point at W 2 = 1.4 GeV 2 which is significantly higher. 



flM TT 




flM L 




where fcj is the c.m.s. momentum in the initial state. Note that the single 
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Fig. 2. Left panel: The integrated cross section Q 6 <7 7 * p—^Op (in units of /iftxGeV 6 ) 
as a function of Q 2 for W = 1.11 GeV (lower curve) and VF = 1.15 GeV (upper curve). 
Right panel: The S-wave (solid) vs. the P-wave (dashed) contribution to the structure 
function Ff (W, Q 2 ) as a function of W 2 for Q 2 = 7.14 GeV 2 . 
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Fig. 3. The structure function Ff(W, Q 2 ) as a function of W 2 scaled by a factor 10 3 
compared to the SLAC E136 data at the average value Q 2 = 7.14 GeV 2 (left panel) and 
Q 2 = 9.43 GeV 2 (right panel). 



Note that in our approximation there is no D-wave contribution, and 
the final state interaction is not included. Both effects can increase the cross 
section so that we consider the agreement as satisfactory. We believe that 
the structure function at W 2 = 1.4 GeV 2 already contains a considerable 
D-wave contribution and also one from the tail of the A-resonance and thus 
cannot be compared with our model, at least in its present form. 

To avoid misunderstanding we stress that the estimates of the cross sec- 
tions presented here are not state-of-the-art and are only meant to provide 
one with the order-of-magnitude estimates of the threshold cross sections 
that are to our opinion most interesting. These estimates can be improved 
in many ways, for example taking into account the energy dependence of 
the generalised form factors generated by the FSI and adding a model for 
the D-wave contributions. The model can also be tuned to reproduce the 
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existing lower Q 2 and/or larger W experimental data. 
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